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We describe the spin dynamics of an arbitrary localized impurity in an insulating 
two dimensional antiferromagnet, across the host transition from a paramagnet 
with a spin gap to a Neel state. The impurity spin susceptibility has a Curie- 
like divergence at the quantum-critical coupling, but with a universal, effective 
spin which is neither an integer nor a half-odd-integer. In the Neel state, the 
transverse impurity susceptibility is a universal number divided by the host spin 
stiffness (which determines the energy cost to slow twists in the orientation of 
the Neel order). These, and numerous other results for the thermodynamics. 
Knight shift, and magnon damping have significant applications to experiments 
on layered transition metal oxides. 
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The recent growth in the study of quasi two dimensional transition metal oxide com- 
pounds (|7|j with a paramagnetic ground state and an energy gap to all excitations with 
a non-zero spin (the 'spin gap' compounds, like SrCu203, CuGeOa and NaV205) has led 
to fundamental advances in our understanding of low dimensional, strongly correlated elec- 
tronic systems. These systems are insulators and are so not as complicated as the cuprate 
high temperature superconductors (which display a plethora of phases with competing mag- 
netic, charge, and superconducting order), and this simplicity has clearly exposed the novel 
characteristics of the collective quantum spin dynamics. 

One of the most elegant probes of these spin gap compounds is their response to inten- 
tional doping by non-magnetic impurities, like Zn or Li, at the location of the magnetic ions. 
Such experiments were initially undertaken on the cuprate superconductors ([|, |^, but their 
analogs in the insulating spin gap compounds have proved to be a most fruitful line of inves- 
tigation They have demonstrated a remarkable property of the paramagnetic ground 
state of the host compound: each non-magnetic impurity has a net magnetic moment of spin 
1/2 located in its vicinity (for the case where the host compound has magnetic ions with spin 
1/2). The confinement of spin is a fundamental defining property of the host paramagnet, 
and is a key characterization of the quantum coherent manner in which the host spins form 
a many-body, spin zero ground state: this confining property was predicted theoretically 
(|^, for the paramagnetic states of a large class of two dimensional antiferromagnets. 

We will describe the quantum theory of an arbitrary localized deformation in such anti- 
ferromagnets; examples of deformations are (i) a. single non-magnetic impurity, along with 
changes in the values of nearby exchange interactions; (ii) change in sign of a localized group 
of exchange interactions from antiferromagnetic to ferromagnetic. Our main concern will be 
the behavior of the impurity as the host antiferromagnet undergoes a bulk quantum phase 
transition from a paramagnet to a magnetically ordered Neel state: we will show that the 
spin dynamics of any deformation is universally determined by a single number — an integer 
or half-odd-integer valued spin S. 

Apart from applications to experiments on materials intentionally driven across a quan- 
tum phase transition, our results also lead to new insights and predictions for the behavior 
of impurities in existing spin gap compounds. The traditional view of the spin gap param- 
agnet is based on strong local singlet formation between nearest-neighbor spins (see Fig |I]A 
below); the resulting picture of doping by a non- magnetic impurity is that the partner spin 
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of tlie impurity site is essentially free. To obtain any non-trivial dynamics one performs 
an expansion about such a decoupled limit, and this yields simple localized spin behavior 
with non-universal details depending upon the specific microscopic couplings. In practice, 
however, spin gap systems are usually well away from the local singlet regime, and strong 
resonance between different singlet pairings leads to appreciable spin correlation lengths: 
their spin gap. A, is significantly smaller than J, a typical nearest-neighbor exchange. A 
systematic and controlled approach for analyzing such a fluctuating singlet state, which we 
advocate here, is to find a quantum critical point to a magnetically ordered state somewhere 
in parameter space, and to then expand away from it into the spin gap state. As we shall 
discuss below, the coupling between the bulk and impurity excitations becomes universal in 
such an expansion, and all dynamical properties depend only upon the bulk parameters, A 
and a velocity c (defined below). 

For clarity, we will state our main results in the context of a simple, explicit theoret- 
ical model; however, they are more general, and apply quantitatively to a broad class of 
experimentally realizable systems. We begin by reviewing the properties of the regular 
antiferromagnet described by the Hamiltonian ^ 

7^ = J 5] S, ■ S, + A J ^ S, ■ S, (1) 

where Sj are spin-1/2 operators on the sites of the coupled-ladder lattice shown in Fig |l|, 
with the A links forming two-leg ladders while the B links couple the ladders. The ground 
state of H depends only on the dimensionless coupling A, and we will restrict our attention 
toJ>0, 0<A<1. AtA = the ladders are decoupled, and each forms a spin singlet 
quantum paramagnet (Fig|l|A). This paramagnetic state continues adiabatically for small 
non-zero A until the quantum critical coupling A = Ac ~ 0.3, where the spin gap vanishes 
as A ~ (Ac — A)"^, where z/ is a known exponent (Qj (the symbol ~ indicates the two 
quantities are asymptotically proportional). For A > Ac, the ground state has long range 
Neel order (Fig|l|B) characterized by the non-zero spin stiffnesses, Psx,psy, which determine 
the energy cost to twists in the order parameter orientation in the x, y directions (we also 
define ps = {psxPsyY^"^)- The low-lying excitations above the Neel state are spin- waves 
which travel with velocities c^, Cy in the x, y directions (with c^/c^ = Psxl Psy\ we define 
c = (ca;Cy)^/^). As A approaches the critical value Ac from above, all the stiffnesses vanish as 
(A — Ac)"^, while the velocities remain finite and non-critical. 
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Introducing a non- magnetic impurity in Ti by removing the spin at site i = X (Fig 
the modified Hamiltonian Tix has the same form as Ti but all links connected to site X do 
not appear in the sums in Eq |I|. The system can be probed by examining its total linear 
susceptibility (x) to a uniform magnetic field, H (under which the Hamiltonian becomes 
Ti-x — gf^B J2i^x H- Sj where fiB is the Bohr magneton and g is the gyromagnetic ratio of the 
ion). This susceptiblity may be written as x = {gfJ'B)^{AXb + Ximp) where A is the total area 
of the antiferromagnet, xt is the bulk response per unit area of the antiferromagnet without 
the impurity, and Ximp is the additional contribution due to the non-magnetic impurity. We 
will now describe the behaviors of Xb and Ximp as the temperature T — > 0, and A moves 
across Ac. 

In the quantum paramagnet, A < Ac, the presence of the spin gap implies that the bulk 
response is exponentially small, Xb = {^/T^Ti^c^)e~^/^^'^ The confinement of a magnetic 
moment in the vicinity of the impurity site implies that there will be Curie like contribution, 
and so 

_ 5(5 + 1) 

where 5 = 1/2 for the model under consideration here (Qj; for a general local deformation, 
we consider Eq ^ as the definition of the value of 5, which, naturally, must be an integer or 
a half-odd-integer. These expressions for Xb and Ximp are exact as T ^ for all < A < Ac. 
Another way of characterizing the confinement of the magnetic moment near X is by looking 
at the time autocorrelation function of a spin at a site i = Y close to X (say, its nearest 
neighbor); at T = this obeys 

Jim (Sy(r) ■ Sy(0)) = ^ 0, (3) 

where r is imaginary time, and my is the local remnant magnetic moment on site Y , which 
is usually significantly smaller than the total impurity moment S appearing in Eq ^ 

Next, we turn to the behavior as T ^ at the critical point A = Ac (more generally, the 
T > results here will apply for A < T < J (p^ < T < J) for A < Ac (A > Ac)). We expect 
that as the spin gap in the quantum paramagnet disappears, the bulk magnon excitations 
will proliferate and their screening will eventually quench the impurity moment — so my 
approaches as A approaches Ac from below. We can anticipate a power-law decay of the 
spin autocorrelations (|^, [71|j , with 

(Sy(r)-Sy(0))~l/r''' (4) 



for large r, T = 0, and A = Ac, and our result for the new, universal exponent rj' is given 
below; standard scaling arguments also imply that my vanishes as my ~ (Ac — X)'^''^^'^. 
The behavior at the critical point therefore appears analogous to that in the overscreened 
multichannel Kondo problem ( |7^ , [7^; in that case, the impurity spin is screened by a bath 
of conduction electrons carrying multiple 'flavors', and also exhibits a power-law decay in its 
autocorrelation. Furthermore, in the multichannel Kondo case, the T dependence of Ximp is 
given essentially by the Fourier transform of Eq ^, that is by Ximp ~ 

This result 

is a consequence of a 'compensation' effect (filj, as the magnetic response of the screening 
cloud of conduction electrons is negligible: the local Fermi levels of up and down electrons 
adjust themselves to the local magnetic field, and hence the susceptibility is not very different 
from the bulk susceptibility except in the immediate vicinity of the impurity spin f|75y. In 
more technical terms, Ximp vanishes in the strict continuum limit, and corrections to scaling 
have to be considered, which lead eventually to Ximp ~ j'-i+V Qur computations show 
that the behavior of Tix at A = Ac is dramatically different: the magnon excitations do not 
have an exact compensation property, and their response is non-zero already in the scaling 
limit. So in a sense, the present problem is simpler than the overscreened Kondo case, and 
naive scaling arguments always work, without inclusion of irrelevant operators — the scaling 
dimension of x is that of inverse energy and so we have one of our central results: 

at A = Ac, where Ci is a universal number independent of microscopic details (as are all the 
Ci introduced below). We computed Ci in the expansion in e = 3 — c/, where d is the spatial 
dimension, and obtained 



_s{s + r] 
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the omitted higher order corrections in Eq ^ will, in general, depend upon S. Comparing 
with Eq ^ we can define an effective impurity spin, Sefr, at the quantum-critical point by 
Ci = 'S'efr(5'efr + l)/3; it is evident that S^q is a universal function of S", is neither an integer 
nor a half-odd- integer, and is almost certainly irrational at e = 1. Also notice that the 
leading corrections in the e-expansion are quite large, and this will be a feature of all the 
results obtained below; accurate numerical estimates require some resummation scheme, but 
we will not discuss this here. For completeness, let us note that at A = Ac, the bulk response 



Xb = C2{kBT)/{hcY, a T-dependence that is also different from the bulk response in 
the overscreened Kondo problem. 

Finally, we describe the situation for A > Ac- Now the presence of Neel order at T = 
implies that the response is anisotropic. Parallel to the Neel order, there is a total magnetic 
moment quantized precisely at S (Qj, and this does not vary under a small longitudinal 
field (there is also a staggered local moment in zero field, as defined by Eq |^, which obeys 
my ~ |A — Xcl^'"^'^). Orthogonal to the Neel order, there is a linear response to a transverse 
field, x±- Foi' the bulk response, we have the well-known result that Xb± is proportional 
to the spin stiffness, Xbi. = Ps/{fi-cY- In contrast, the same scaling arguments leading to 
Eq ^ imply that Ximp± is inversely proportional to Ps, the latter being the only energy scale 
characterizing the ground state as A approaches Ac from above; so another key result of this 
paper is 

Ximp± • (7) 

Ps 

In general d, this relationship is Ximp± = C^ihcY'^^'^^^^'^^^y pl^^'^^^\ and the e-expansion of C3 
is 

155 / 1 15^+1 



, , 1 - (1.193 + 0.5535 + 0.4195^)e + C(e^) , (8) 

v22 V 2e / L J 

where Sd = 2/ (r((i/2)(47r)'^/^). Note that ps vanishes, and so Ximp± diverges, as A approaches 

Ac. Turning to T > but very small, in = 2 and in the absence of any spin anisotropy, 

strong angular fluctuations cause the Neel order to vanish at any non-zero T. Then the 

susceptibility takes the rotationally averaged value Ximp = 

Sy{3kBT) + (2/3)ximp±, where 
the first term is the contribution of the net moment noted earlier (note that this term has a 
coefficient 5^ and not S{S + 1), because the locking of the moment orientation to the local 
Neel order makes it behave classically). In practice, this averaged Ximp will not be observable 
as even an extremely small anisotropy will pin the Neel order below a small T > 0. Our 
results for x are summarized in Fig |^. 

The next two paragraphs contain a technical interlude which outlines the field-theoretic 
derivation of the results above — details appear elsewhere (fT^j. We describe the bulk-ordering 
transition by a (i + 1-dimensional field theory with action Si, of a field (pai^, t) (a = 1 ... 3) 
representing the collinear Neel order parameter (Qj. This is coupled by the action iSimp to an 
impurity spin at a; = with orientation given by the unit vector ria- The partition function 
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is / V(j){x, T)Vn{T) exp(— iSf, — iSimp) with 



Sh = Id xdr 



<Simp — I dr 



iSAain)^^ - '-foSna{T)(t)a{x = 0, r) 
ar 



(9) 
(10) 



where eap-yOAp/dn.^ = na, and the term proportional to A{n) is a Wess-Zumino form repre- 
senting the Berry phase of the impurity spin. The bulk transition in Sb is driven by tuning 
the coupling r through a critical value r^, which therefore plays a role similar to A; the A < Ac 
(A > Ac) region of the lattice antiferromagnet Ti. maps onto the r > Tc (r < Tc) region of the 
field theory Sb- Quite generally, any local deformation of the antiferromagnet is described 
by the action Sb + Si^p, where S, defined as the integer or half-odd- integer appearing in 
Eq ^, is (roughly) the net local imbalance of spin between the two sublattices. Changes 
in exchange constants lead to additional terms like / (ir0^(x = 0,r) which are all strongly 
irrelevant under the renormalization group (RG) analysis in powers of e. The r = 0, Qq = 
case of Eqs Pi|To| was considered earlier by Sengupta ( [7^) (and related models in |II|)) in a 
non-local formulation in which (f)a{x 7^ 0, r) was integrated out: however, such a model has a 
pathological response to even an infinitesimal field H (the energy is unbounded below), and 
the quartic go coupling is essential to stabilize the system, and to all the results obtained 
here. Further, the local formulation here facilitates development of the RG to all orders. 

The RG analysis of Sb + Si^p is carried out by the methods of 'boundary critical phe- 
nomena' of a ((i + l)-dimensional system with a 1-dimensional 'boundary' at x = 0, 
which constitutes a 'dimensional reduction' oi d > 1 (contrast this with the case of a. {d+ 1)- 
dimensional system with a d-dimensional boundary, with a dimensional reduction of 1, 



which has been invariably ( fZ^ [7^ j considered earlier, as in all the Kondo problems). The 
irrelevance of the boundary 'mass' term 0^(x = 0, r) implies that there is only an 'ordi- 
nary' transition at the position of the bulk critical point (this has been implicit in 
our earlier discussion), and there are no analogs of the 'surface', 'special', and 'extraor- 
dinary' transitions j^TSj). The RG analysis of the bulk action Sb is now standard text- 
book material — we will not reproduce it here, and will follow the notation of d^ ). We 
introduce renormalized fields = V^0k, n = \fZ'n^^ and renormalized couplings by 



fi'o = (/^Vc)(^4/^^'S'd+i)(7, 7o = {^"cY^'^^Z^i/ \l ZZ'Sd+i)'^ where /i is a renormalization in- 
verse length scale, Sd = T{d/2 — l)/(47r'^/^), and the bulk renormalization factors Z, are 



specified in ^^). For the new boundary renormalization factors, we obtained to two loops 



Z' = l-27Ve + 7Ve and Z., = l + 7r^{S{S +1) -l/3)g-f^ /{6e). These lead to the (3 function 



for g found in and the new f3 function for the boundary couphng 

Pil) = + 7' - 7' + ^ + tis{s + 1) - 1/3)^77=^ + o ((7, VaY) (11) 



2 ' ' 144 3 

The critical fluctuations at the boundary are therefore controlled by the flxed point values 
1 = 1*1 9 = 9* (both nonzero) at which both j3 functions vanish, and canonical methods 
then imply the exponent 
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— + — (S(S + l)-l/3) 
242 11 ^ ^ ' ' ' 



t^ + 0{€'). (12) 



Eq ^ can now be obtained by the methods of (]^^, while Eq ^ follows directly from a 
renormalized perturbation theory in the ordered phase at T = 0. We conclude our technical 
interlude by noting that our RG scheme shows directly that the only graphs which contribute 
to the renormalization of 70, beyond those arising from the wavefunction renormalization 
Z\ must include a factor of the bulk interaction qq] this implies that Z^ = 1 for (7 = 0, and 
shows that for the models of (j^ [J^ j the one-loop exponent r^' = e is exact. 

The above methods can be extended to determine the behavior of other observables in 
the regimes of Fig |^. We mention a few: 

(z) Entropy: In the paramagnetic phase (A < Ac) there is clearly a residual entropy of 
ln(2S'+ 1) as T — > 0. At A = Ac, the e-expansion shows that this is modifled to ln(2S' + 1) — 
S{S + l)(33e/160)i/2 _^ 0(^3/2)^ in the Neel state (A > Ac, the Neel order pinned by 

some small spin anisotropy) the impurity entropy vanishes as T'^ at low T. 
{a) Knight shift: We restrict the discussion here to the intermediate quantum-critical region 
of Fig^, T > |A — Acl'^. The shift in the NMR resonance frequency is proportional to the local 
response in the presence of a uniform external fleld, x{^)- the vicinity of the impurity 
{e.g. at site i = Y) x{^) ~ Well away from the impurity (|x| 00), apart from 

the bulk response of the antiferromagnet, there are staggered and uniform contributions 
which decay exponentially with |x| on a scale ~ he/ {^kBT). 

[iii) Magnon damping: In the quantum paramagnet (A < Ac), and at T = 0, the pure 
antiferromagnet has a pole in the dynamic spin structure factor ~ 1/(A — hu) at the 
antiferromagnetic ordering wavevector from the triplet magnon excitations. In the presence 
of a dilute concentration of impurities, rii, this pole will be broadened on an energy scale 
F; scaling arguments and the structure of the flxed point found here imply the exact form 
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(jl^ r ~ ni{hcYA^~'^. We argue that this damping mechanism is the main ingredient in the 
broadening of the 'resonance peak' observed recently in Zn-doped YBa2Cu307 (^4\)- Using 
the values he = 0.2a eV (a is the lattice spacing), A = 40 meV, and = 0.005/a^, we 
obtain the estimate F = 5 meV, which is in excellent accord with the observed linewidth of 
4.25 meV We have also studied the lineshape of the magnon peak (fT^, and find that 
it is asymmetric at very low T, with a tail at high frequencies: it would be interesting to 
test this in future experiments. 

We have described the highly non-trivial, collective, quantum spin dynamics of a single 
impurity in a strongly correlated, low dimensional electronic system. The problem maps onto 
a new boundary quantum field theory, Eqs p|,p!0|, and is therefore also of intrinsic theoretical 
interest: unlike previously studied quantum impurity problems, there is a complicated inter- 
ference between bulk and boundary interactions, and its proper description is the key to the 
physical results we have obtained. Our theoretical results for the magnon damping in the 
spin gap phase are in good agreement with existing experiments ^^). Studies of materials 
exhibiting other aspects of the regimes of Fig ^ appear possible, and we hope they will be 
undertaken; spin gap compounds can be driven across the transition by, say, application of 
hydrostatic pressure, or by doping with other impurities which have the same spin as the 
host ion they replace and do not change the sign of the exchange constants Quantum 
Monte Carlo simulations should also allow more accurate determination of the universal 
constants Ci and C3. 
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FIG. 1: The coupled ladder antiferromagnet. The A links are full lines and have exchange J, while 
the B links are dashed lines and have exchange AJ. The paramagnetic ground state for A < Ac 
is schematically indicated in (A): the ellipse represents a sing let valence bond, (| Ti> - I iT))/\/2 
between the spins on the sites. The Neel ground state for A > Ac appears in (B). 
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FIG. 2: The impurity Hamiltonian 7ix in which the spin and hnks on site i = X have been 
removed. 
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FIG. 3: Summary of the results for the bulk and impurity susceptibilities of TCx- The constants 
are universal numbers, insensitive to microscopic details like variations in the magnitude or sign of 
the exchange constants in the vicinity of the impurity, or presence of additional, nearby, vacancies 
or impurity ions with different spins. The constants Ci and C3 depend only on the integer /half- 
odd-integer valued S, and we can view Eq |2|, the T ^ limit of Ximp in the paramagnet (A < Ac), 
as the experimental definition of S. For the case in which non-magnetic impurities are added in a 
localized region, with no modification of exchange constants, S is the net imbalance of spin between 
the two sublattices. The constant Ci defines the effective spin at the quantum-critical point by 
Ci = Ses{SeS + l)/3, and Ses is neither an integer nor a half-odd-integer. 
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